Dynamics of Vacillating Voters 
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We introduce the vacillating voter model in which each voter consults two neighbors to decide its 
state, and changes opinion if it disagrees with either neighbor. This irresolution leads to a global 
bias toward zero magnetization. In spatial dimension d > 1, anti-coarsening arises in which the 
linear dimension L of minority domains grows as t^/f'^+i). One consequence is that the time to 
reach consensus scales exponentially with the number of voters. 
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The voter model [ij gives an appealing, albeit ideal- 
ized, description for the opinion dynamics of a socially 
interacting population. In this model, each node of a 
graph is occupied by a voter that has one of two opin- 
ions, f or j. The population evolves by: (i) picking a 
random voter; (ii) the selected voter adopts the state of a 
randomly-chosen neighbor; (iii) repeating these steps ad 
infinitum or until a finite system necessarily reaches con- 
sensus. Descriptively, each voter has no self confidence 
and follows one of its neighbors. With this dynamics, 
a voter chooses a state with a probability equal to the 
fraction of neighbors in that state, a feature that renders 
the voter model soluble in all dimensions P, [1] . 

In this work, we investigate a variation that we term 
the vacillating voter model. By vacillating, we mean that 
a voter very much lacks confidence in its state. In an 
update, if a voter happens to select a random neighbor 
of the same persuasion, the voter is still not convinced 
that this state is right. Thus the voter selects another 
random neighbor and adopts this state. This vacillation 
causes a voter to change state with a larger probability 
than the fraction of disagreeing neighbors, and leads to a 
bias toward the zero-magnetization state in which there 
are equal densities of voters of each type. 
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FIG. 1: Illustration of an update for the vacillating voter on 
the square lattice (left and middle). For the configuration on 
the right, the central voter flips with probability 5/6 because 
out of the 6 ways of selecting two neighbors, only one choice 
leads to both neighbors agreeable (dashed). 

Thus vacillation inhibits consensus, but due to a dif- 
ferent mechanism than that in the prototypical Axelrod 
model (sl, the bounded compromise model [1] and its 
variants [5I|. For these latter models, consensus is hin- 



dered because of the absence of interaction whenever two 
agents become sufficiently incompatible. For vacillating 
voters, it is individual uncertainty that forestalls consen- 
sus. The vacillating voter model also differs from models 
that incorporate "contrarians" [6|] because voters still try 
to imitate their neighbors. 

The update steps in the vacillating voter model are: 

1. Pick a random voter. 

2. The voter picks a random neighbor. If the neighbor 
disagrees with the voter, the voter changes state. 

3. If the neighbor and the voter agree, the voter picks 
another random neighbor and adopts its state. 

4. Repeat steps 1 and 2 ad infinitum or until consensus 
is reached. 

For example, the probability that a vacillating voter on 
the square lattice fiips is 0, ^, |, and 1, respectively, when 
the number of anti-aligned neighbors is 0, 1, 2, and > 3 
(Fig. [T]). In contrast, for the classic voter model, the fiip 
probability is |, where k is the number of neighbors of 
the opposite opinion. We now explore the consequences 
of this vacillation on voter dynamics. 

Consider first the mean-field limit. Here the density x 
of t voters obeys the rate equation 
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The first term on the right accounts for the loss of t voters 
in which a t voter is first picked (factor x), and then 
the neighborhood cannot consist of two t voters (factor 
1 — x^). Similarly, in the second (gain) term, a J, voter is 
first picked, and then the neighborhood must contain at 
least one t voter. The factorized form shows that there 
are unstable fixed points at a; = 0, 1 and a stable fixed 
point at X = 1/2. Thus a population is driven to the 
zero-magnetization state. 

However, because consensus is the only absorbing state 
of the stochastic dynamics, a finite population ultimately 
reaches consensus. To characterize the evolution to this 
state, we first study the exit probability f„, defined as 
the probability that a population of N voters ultimately 
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reaches t consensus when there are initially n t voters. 
Then obeys the backward equation Q 



n+1 £n+l 



£n, (2) 



where Wn^m is the probability for the transition from 
the state with n t voters to to "f voters in an update. 
This equation expresses the probability to exit from n 
as the probability to take one step (the factors w) times 
the probability to exit from the point reached after one 
step. In the large-7i limit, we write x — n/N, and the 
transition probabilities become 

Wn^n+l = (1 - X) [l - (1 - X)^] 
Wn^n = X^ + (1 - X)^. 

Substituting these in writing £n±i Six ± Sx), 
and expanding to second order in 6x, gives 

3x{l-x) d^£ 

w 

with solution 
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£{x) = 



x-l/2 



1/2 




(4) 



Notice that S{x) approaches the constant value 1/2 for 
increasing N (Fig. (2]), reflecting the bias towards the zero- 
magnetization state. Almost all initial states are driven 
to the potential well at x = 1/2, so that the exit prob- 
ability becomes independent of the initial density of | 
voters. 




FIG. 2: Exit probability £{x) versus the density of | voters x 
for the case = 16, iV = 25 and N = 100. 



Similarly, we study the time to reach consensus as a 
function of the initial composition of voters. Let i„ de- 
note the time to reach consensus (either all t or all |) 
when starting with n | voters in a population of TV vot- 
ers. Similar to tn obeys the backward equation Q 
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where 6t = 1/N is the time elapsed in an update. In the 
large-n limit, this equation becomes 
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The formal solution is again elementary, but the result 
can no longer be expressed in closed form. The main 
result is that the consensus time scales as e°^, with a 
a constant of order 1. In contrast to the classical voter 
model, the global bias drives the system into a potential 
well that must be surmounted to reach consensus. Thus 
the consensus time is anomalously long. 

In one dimension, a voter changes its opinion if at least 
one of its neighbors is in disagreement. For example, a | 
voter flips with rate 1 if the neighborhood configurations 
are ttJ-j itt: and |ti- As an amusing side-note, this dy- 
namics is equivalent to rule 178 of the one-dimensional 
cellular automaton Q, except that this rule is imple- 
mented asynchronously in the vacillating voter model. 
In the framework of the Ising-Glaubcr model @, the flip 
rate of a voter at site i, whose states are now represented 
by a., = ±1, is 



[o'i(fi+l+crj-l) + o'i_icri+i —3] 



(7) 



with {cr} denoting the state of all voters and {a'}i the 
state where the i"^ voter flips. The first two terms corre- 
spond to conventional Glauber kinetics, but as mentioned 
parenthetically in Ref. 9], the presence of the ai-iUi+i 
term couples the rate equation for the mean spin to 3- 
body terms and the model is not exactly soluble. 

The mean spin, Sj = {aj) = Xl{cr} '^jPii'^}'^^) evolves 
according to 

^ = 5]a, [^^{a'}.^ W)F({a'}.;t) 



i'y} 



,{{a}^{a'},)P{{a}-t) , (8) 



which reduces to, after straightforward but tedious steps, 
ds 1 

= 2 ^^'^^ ^ ^'^^ ^ {<^3-i<^3'^i+i) - 3sj) . (9) 

In a similar spirit, the rate equation for the nearest- 
neighbor correlation function, (cTjaj+i), is 

^^^^^ = \ [{'y3-i{<yo+<^j+i)) + {{<^'j+'^m)'^'m)\ 

+ l-3(a,a,+i) (10) 

We can simplify Eq. (jlOp by considering domain 
walls — nearest-neighbor anti-aligned voters — whose den- 
sity is given by p = (1 — {(liai+i)) /2. According to the 
flip rate in Eq. ([7]) , an isolated domain wall diffuses freely, 
just as in the pure voter model. However, when two do- 
main walls are adjacent, they annihilate with probability 
1 /3 or one hops away from the other with probability 2 /3. 
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This process is isomorphic to single-species annihilation, 
A 0, but with a reduced reaction rate compared to 
freely diffusing reactants because of the nearest-neighbor 
repulsion. The domain wall density still asymptotically 
decays as t^^^^ with an amplitude that depends on the 
magnitude of the repulsion. 

Because domain walls are widely separated at long 
times, the second-neighbor correlation function is 

(CTjaj+2) = +prob(0 or 2 walls between j and j + 2) 
— prob(l wall between j and j + 2) 
w 1 - 2p. 

Using the approximation of widely separated domain 
walls, {a j a ~ i'^j'^j+i) = ™2, and the rate equa- 
tion for nearest-neighbor correlation function m2 be- 
comes = 1 — 7712 , with solution 



m2(i) = 1 -I- [777(0)^ - 1] e+ 



(11) 



Here we chose the uncorrelated initial condition, so that 
7772(0) = 777(0)^, where 7^(0) = {sj{0)) is the average 
magnetization at t — 0. 



Since the density of f voters is a; = (1-1- 77i)/2, while 
771(00) = 2£{x) — 1, the exit probability £{x) becomes 



£{x) = - [(22; - l)e2^(^-i) + 1 



(15) 



This result is in excellent agreement with our simula- 
tion results (Fig. [3|). For small systems {N — 25 and 
36), we directly measure the probability £{n) that the 
population ultimately reaches a t consensus when there 
arc initially n | voters and averaged over 5000 realiza- 
tions of the dynamics. We also verified Eq. for large 
systems {N — 1000 nodes) by a different approach that 
avoids the need to measure £{n) directly by simulating 
until ultimate consensus. Instead, we run the dynamics 
up to 1000 time steps and measure the magnetization 
at this time. We then average over 200 realizations of 
the process to obtain 777(00) and finally obtain £{x) from 
£{x) = (l-|-77i(oo))/2. We again find excellent agreement 
with our prediction (jl5p . 
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FIG. 3: Exit probability £{x) as a function of the initial den- 
sity of I voters x for a one dimensional system composed of 
25, 36 and 1000 voters respectively The voter model result, 
£{x) = X, that follows from magnetization conservation is 
shown for comparison. 

Let us now return to the rate equation ^ for the mean 
spin. For a spatially homogeneous system, (sj) are all 
identical and the magnetization is 777 = (sj). Also, we 
follow Ref. [13] and decouple the 3-spin correlation func- 
tion as (aj-iajaj+i) « 7777772. Then by averaging over all 
sites, the rate equation equation ([9]) becomes 

dm 1 



-(777 7r72 — 771) 



(777(0)^-1), (12) 



dt 2' ~^ 2 
whose solution, for the initial condition 777(0), is 

777(t) - 777(0) e^(l-^")(™(0)'-l) . 



(13) 



Thus we obtain a non-trivial relation between final mag- 
netization 777(00) and 77i(0) 

777(0) e5("(0)'-i) 



7^7(00) 



(14) 
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FIG. 4: Double logarithmic plot of the number of | voters 
versus time on the square lattice starting from a 4 x 4 square 
of I voters in a background of J, voters. 

The vacillating voter model in greater than one dimen- 
sion has the new qualitative feature that small minority 
domains tend to grow. This anti-coarsening is a manifes- 
tation of the bias toward the zero-magnetization state. 
To appreciate how this anti-coarsening arises, consider a 
circular two-dimensional island domain of f voters of lin- 
ear dimension L and area A in a sea of J, voters. For large 
L, each voter at the interface has the same local environ- 
ment, so that there is no environmental bias. However, 
there are slightly more | voters just outside the circle 
that t voters just inside. In a time of the order oi 6t ^ L 
each interface voter is updated once, on average, so that 
the island area increases by an amount SA that is of the 
order of the difference in the number of f and J, voters at 
the interface. Thus ^ ~ -^j which gives L ^ t^/'^. In d 
dimensions, this same reasoning gives 
probed for this anti-coarsening by simulating the evolu- 
tion of an initial small square domain of "f voters in a | 
background in two dimensions (Fig. 3]). Although such 
domains do not remain contiguous, the data suggest that 
the number, or occupied area, of f voters grows as t"', 
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with a around 0.73, in reasonable agreement with our 
expectation a = 2/3. 




FIG. 5; Exit probability £(x) as a function of the initial den- 
sity of I voters x for a square lattice of 16, 25, 36 and 49 
voters, respectively, with periodic boundary conditions. 





the system achieves consensus without first being drawn 
to this attractor. Thus the exit probability S{x) should 
be nearly independent of x for almost all a;, just as in 
the mean- field limit. Simulations of the vacillating voter 
model on the square lattice (Fig. [5]) confirm that £{x) 
approaches 1/2 for a progressively wider range of x as 
L increases. Simulations also show that the correlation 
function Ci = {ai^jai^j+i) does not approach 1 in the 
long-time limit, as in one dimension or in the pure voter 
model in two dimensions. Rather, Ci reaches the sta- 
tionary value 0.31, so that domains of opposite opinions 
coexist (Fig. [S]), and only a rare macroscopic fluctuation 
allows consensus to be reached. 

In summary, when vacillation is incorporated into the 
voter model, consensus is inhibited but not prevented. In 
the mean-field limit, the vacillation drives a population 
away from consensus and toward the zero-magnetization 
state. A finite system ultimately achieves consensus only 
via a macroscopic fluctuation that allows the system to 
escape this bias-induced potential well. Because of the 
bias, the probability to reach f consensus is essentially 
independent of the initial composition of the popula- 
tion. In one dimension, the system coarsens, albeit more 
slowly than in the pure voter model because of the re- 
pulsion of neighboring domain walls, and the probability 
to reach the final state of t consensus has a non-trivial 
initial state dependence. In two and higher dimensions, 
domains slowly anti-coarsen to drive the system to the 
zero-magnetization state. The overall behavior is quali- 
tatively similar to that of the mean-field vacillating voter 
model, and very different from the pure voter model. 
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FIG. 6: Snapshots of the vacillating (left) and pure (right) 
voter model on a 50 x 50 lattice starting with a random zero- 
magnetization state after 100 time steps. The correlation 
function Ci equals 0.31 (left) and 0.59 (right) respectively. 

A system with non-zero initial magnetization is there- 
fore again drawn to the attractor where the density x of 
t voters equals 1/2 before final consensus is eventually 
reached. It is only for x initially very close to or 1 that 
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